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ABSTRACT 

We present the final installment of an intensive 13-year study of variations of the 
optical continuum and broad H/3 emission line in the Seyfert 1 galaxy NGC 5548. 
The data base consists of 1530 optical continuum measurements and 1248 H/3 
measurements. The H/3 variations follow the continuum variations closely, with a 
typical time delay of about 20 days. However, a year-by-year analysis shows that the 
magnitude of emission-line time delay is correlated with the mean continuum flux. We 
argue that the data are consistent with the simple model prediction between the size of 
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the broad-line region and the ionizing luminosity, r oc L^ n . Moreover, the apparently 
linear nature of the correlation between the H/3 response time and the nonstellar 
optical continuum F opt arises as a consequence of the changing shape of the continuum 
as it varies, specifically F opt cx F^y 6 . 

Subject headings: galaxies: active — galaxies: individual (NGC 5548) — galaxies: 
nuclei — galaxies: Seyfert 

1. Introduction 

The nature of the broad-line region (BLR) is arguably one of the major remaining mysteries 
in active galactic nuclei (AGNs). Whereas the supermassive black hole/accretion-disk paradigm 
has become increasingly more secure, there is still no consensus about the origin of the broad 
emission lines that are prominent features of the UV/optical spectra of these sources. This is 
not to say that progress has not been made. Years of spectroscopic study have uncovered a rich 
phenomenology (see Sulentic, Marziani, & Dultzin-Hacyan 2000 for a recent review). 

Of particular importance has been the recognition that the emission-line fluxes vary in 
response to continuum variations, with a small time delay (days to weeks for Seyfert galaxies) due 
to light travel-time effects within the BLR. Well before these time delays were first accurately 
measured, this led to the seminal paper on "reverberation mapping" (Blandford & McKee 1982), 
a tomographic method of determining the structure and kinematics of the BLR. While the full 
potential of reverberation mapping has yet to be realized, use of this technique has provided BLR 
sizes in approximately three dozen AGNs (see compilations by Wandel, Peterson, & Malkan 1999 
and Kaspi et al. 2000). Even more importantly, measurement of the emission- line time delays (or 
"lags"), combined with measurements of the line width, have led to estimates of the masses of the 
central black holes. 

In late 1988, we began a program of spectroscopic monitoring of optical variations in the 
Seyfert 1 galaxy NGC 5548. This activity was organized by an informal consortium called the 
"International AGN Watch" (Alloin et al. 1994; Peterson 1999). Initially, the optical monitoring 
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program was undertaken in support of an ultraviolet monitoring campaign carried out with the 
International Ultraviolet Explorer (IUE; Clavel et al. 1991 [the first paper in the same series as 
this work, and hereafter referred to as Paper I]). The resulting optical data set (Peterson et al. 
1991, hereafter Paper II) proved to be so rich that it was decided to continue this program. Since 
that time, a number of additional contributions to this series (Peterson et al. 1992 [Paper III]; 
Peterson et al. 1994 [Paper VII]; Korista et al. 1995 [Paper VIII]; and Peterson et al. 1999 [Paper 
XV]) have described our continuing optical monitoring program. A complete summary of the work 
on NGC 5548 by this consortium appears in Table 1. Reviews of the progress of this program are 
also available (e.g., Peterson 1993, 2001; Netzer & Peterson 1997). 

This paper represents the final installment of our 13-year program of optical monitoring of 
NGC 5548. As is customary for the papers in this series, we will focus primarily on the data and 
the salient results of cross-correlation analysis of the continuum and emission-line light curves. 
Further in-depth analysis will be left to subsequent papers. In § 2, we describe the observations, 
data reduction, and intercalibration procedures that we have used to construct a homogeneous 
data base of optical continuum and H/3 emission-line fluxes. In § 3, we describe the time-series 
analysis that we have undertaken to determine the time scale for response of H/3 to continuum 
variations. In § 4, we discuss some of the implications of our study. Our results are summarized in 
§ 5. We also note that some of the data presented here have already appeared in other contexts, 
specifically as part of a short but intensive monitoring program carried out in 1998 June (Dietrich 
et al. 2001) and as part of a long-term photometric monitoring campaign (Doroshenko et al. 2001). 

2. Observations and Data Analysis 

2.1. Spectroscopic Observations 

Here we analyze spectroscopic observations made between 1995 November 22 (Julian Date 
= JD2450044) and 2001 September 21 (JD2452174); UT dates are used throughout this paper. 
During this six-year period, a total of 543 individual spectra were obtained, as summarized in 
Table 2. Column (1) gives a code for each data set; these are the same codes that have been used 
throughout the NGC 5548 monitoring project. A more complete log of observations, as has been 
published in our previous papers on NGC 5548, can be found at the International AGN Watch 
site on the World-Wide WebQ. 

The data obtained during Year 8 of this program (see Table 1) were previously published in 
Paper XV. Here we completely reanalyze the data from Year 8 because 46 additional spectra from 



lr The light curves and complete logs of observation are available in tabular form at URL [ittp://www.astronomy.ohio^ 
5tate.edu/~agnwatch/. All publicly available International AGN Watch data can be accessed at this site, which also 
includes complete references to published AGN Watch papers. 
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set "F" have become available to us since Paper XV was published^. 

The spectroscopic images were processed by individual observers in standard fashion for CCD 
frames, including bias subtraction, flat-field correction, wavelength calibration, and flux calibration 
based on standard-star observations. The flux calibration is then refined by scaling each spectrum 
to a constant value of the [O in] A5007 narrow-line flux, as we have done in each previous paper in 
this series. All spectra are gray-scaled in flux (i.e., adjusted by a constant multiplicative factor) to 
a constant flux of F([0 in] A5007) = 5.58 x 10 -13 ergs s _1 cm~ 2 . This absolute flux was determined 
from spectra made under photometric conditions during the first year of this program (Paper II). 
Scaling is accomplished automatically by use of the van Groningen &: Wanders (1992) spectral 
scaling software, as described in Paper XV. 

After flux scaling, measurements of each of the spectra are made. The continuum flux at 
~ 5100 A (in the rest frame of NGC 5548, z = 0.0167) is determined by averaging the flux in 
the 5185-5195 A bandpass (in the observed frame). The Fl/3 emission-line flux is measured by 
assuming a linear underlying continuum between ~ 4790 A and ~ 5170 A, and integrating the 
flux above this continuum between 4795 A and 5018 A (all wavelengths in the observed frame). 
The long-wavelength cutoff of this integration band misses some of the H/3 flux underneath 
[O in] A4959, but avoids the need to estimate the Fe n contribution to this feature and still gives a 
good representation of the H/3 variability. We also note that no attempt has been made to correct 
for contamination of the line measurement by the narrow-line component of H/3, which is of course 
expected to be constant. 

Even after scaling all of the spectra to a common value of the [O in] A5007 flux, there are 
systematic differences between the light curves produced from data obtained at different telescopes. 
As in our previous papers, we correct for the small offsets between the light curves from different 
sources in a simple, but effective, fashion. We attribute these small relative offsets to aperture 
effects (Peterson et al. 1995), although the procedure we use also corrects for other unidentified 
systematic differences between data sets. We define a point-source correction factor <p by the 
equation 

F(H/3) true = ^F(H/3) observed . (1) 

This factor accounts for the fact that different apertures result in different amounts of light loss for 
the point-spread function (which describes the surface-brightness distribution of both the broad 
lines and the AGN continuum source) and the partially extended narrow-line region. 

After correcting for aperture effects on the point-spread function to narrow- line ratio, we also 
correct for the different amounts of starlight admitted by different apertures. An extended source 



2 We also note in passing that we have included in the AGN Watch archive five additional set "H" spectra from 
Years 3 and 4 (details are provided on the AGN Watch website), although these data have not been incorporated 
into the light curves that are currently available. 



- 5 - 



correction G is thus defined as 

F x (5100 A) true = <^F\(5100 A) 

observed 

G. (2) 

Intercalibration of the individual data sets is then accomplished by comparing pairs of nearly 
simultaneous observations from different data sets to determine for each data set the values of the 
constants ip and G that are needed to adjust the emission-line and continuum fluxes to a common 
scale. Furthermore, the formal uncertainties in ip and G reflect the uncertainties in the individual 
data sets, so we can determine the nominal uncertainties for each data set if we assume that the 
errors add in quadrature. 

As in our previous work, the data are adjusted relative to data set "A" because these 
data are fairly extensive, overlap well with most of the other data sets, and were obtained 
through a reasonably large aperture (5".0 x 7".5). As in Paper XV, fractional uncertainties of 
Ccont/F\(5100 A) rs 0.020 and o"ii nc /F(H/3) 0.020 for the continuum and H/3 line, respectively, 
are adopted for the similar, large-aperture, high-quality data sets "A" and "H," based on the 
differences between closely spaced observations within these sets (see also Peterson et al. 1998a). 
For the other data sets, it was possible to estimate the mean uncertainties in the measurements 
by comparing them to measurements from other sets for which the uncertainties are known and 
by assuming that the uncertainties for each set add in quadrature. 

The intercalibration constants we use for each data set are given in Table 3, and these 
constants are used with equations (||) and @) to adjust the spectral measurements. We note that 
there were a few cases in which the intercalibration constants had to be evaluated on a year-to-year 
basis to effect an acceptable intercalibration with the rest of the data. 



2.2. Photometric Observations 

In addition to the spectroscopic observations, numerous y-band photometric observations of 
NGC 5548 were made. The sources of these observations are summarized in Table 4 and a more 
complete log of observations is available on the AGN Watch website. Many of the observations 
used in this analysis were also included in a study by Doroshenko et al. (2001). 

In order to effect an intercalibration with the spectroscopic light curve, we compared closely 
spaced pairs of spectroscopic and photometric observations and performed a least-squares fit to 
the equation 

F A (5100A) = a + bF v , (3) 

where the magnitude V oc —2.5 log Fy. This fit was based on pairs of observations obtained within 
one day of each other. The raw T/-band fluxes were adjusted using eq. (D to the F A (5100 A) scale 
defined by the spectrophotometry. 



-6- 



2.3. The Merged Light Curves 

A complete continuum light curve was assembled by merging the spectrophotometric and 
F-band photometric light curves described above. The final light curves for the optical continuum 
and H/3 fluxes were then constructed by computing a weighted average of points separated in 
time by less than 0.1 day. The final light curves for Years 8 through 13 are given in Table 5. It 
is important to note that all measurements are in the observer's frame and are uncorrected for 
Galactic extinction. Also, the continuum flux includes a contribution from the host galaxy, which 
we estimate to be 3.4 x 10 _15 ergs s -1 cm -2 A — 1 (Romanishin et al. 1995), and the emission-line 
flux includes the narrow-line component, which we estimate to be ~ 8.2 x 10~ 14 ergs s _1 cm~ 2 
(Wanders & Peterson 1996); both of these values are in the observer's reference frame. 

3. Variability Analysis 

The data given in Table 5 can be combined with our previously published light curves 
(Papers II, III, VII, VIII, and XV), yielding homogeneous light curves that cover a span of 4756 
days. The combined data are shown in Figure 1. These light curves are comprised of 1530 
continuum measurements and 1248 line measurements. In this section, we will summarize the 
basic characteristics of the thirteen-year homogeneous data base. 

3.1. Characteristics of the Data Base 

In Tables 6 and 7, we provide a summary of the basic characteristics of the continuum and 
H/3 emission-line light curves, respectively, as shown in Figure 1. Column (1) indicates the subset 
considered, with time ranges as defined in Table 1. The number of observations in each subset is 
given in column (2), and columns (3) and (4) give the average and median intervals, respectively, 
between successive observations. The mean and root-mean-square (rms) fluxes, F A (5100A) m 
Table 6 and F(H(3) in Table 7, are shown in column (5). Two standard measures of variability, F v3iT 
and -Rmax, are given in columns (6) and (7), respectively. The parameter F var is the rms fractional 
variability, corrected for measurement error, as defined by Rodrfguez-Pascual et al. (1997), and 
-Rmax is simply the ratio of maximum to minimum flux. Both of these parameters are affected by 
contamination of the measured quantities by constant-flux components, namely the underlying 
host galaxy in the case of the continuum, and the narrow H/3 emission component in the case of 
the line. If we use the values for the host- galaxy contribution and the narrow-line H/3 flux given 
in the last section, the values of F var increase to 0.380 and 0.242 for the continuum and H/3 line, 
respectively, and the values of -R max increase to 9.44 ± 1.41 and 6.26 ± 0.62, respectively. It is 
notable that the optical nonstellar continuum in NGC 5548 has varied by an order of magnitude 
over this 13-year monitoring program. 
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3.2. Time-Series Analysis 

As in Paper XV and previous papers, we compute the time delay between continuum variations 
and the H/3 response by cross-correlation of the continuum and emission-line light curves. We 
perform the cross-correlation analysis in two ways, using the interpolation cross-correlation 
function (ICCF) method of Gaskell & Sparke (1986) and Gaskell & Peterson (1987), and the 
discrete correlation function (DCF) method of Edelson & Krolik (1988). In both cases, we use the 
implementation described by White & Peterson (1994). 

The results of the cross-correlation analysis are shown in Figure 2 and summarized in Table 
8. The entries in Table 8 give the subset used in the analysis (column 1), the centroid T ccn t of the 
ICCF (column 2), the peak r pca k of the ICCF (column 3), and r max , the value of the ICCF at 
Tpeak- The value of r ccnt is computed using only those points near the principal peak with values 
exceeding 0.8r max . The uncertainties quoted for r ccnt and r pe ak have been computed using the 
model-independent FR/RSS Monte-Carlo method described by Peterson et al. (1998b). 



4. Discussion 

It was noted in Paper XV that statistically significant year-to-year changes in the H/3 
lag have occurred, and that these changes are correlated with the mean continuum flux. The 
additional data described here confirm this result and allow us to investigate it further. In 
Figure 4, we plot r cent as a function of the mean starlight-corrected continuum flux for each 
year F opt = (Fa(5100A)) — F ga i, where F ga i is our best estimate of the host-galaxy starlight 
contribution through our standard aperture, 3.4 x 10 -15 ergs s _1 cm~ 2 A -1 at 5100 A (Romanishin 
et al. 1995); F opt is thus the mean nonstellar flux from the AGN component only. We add to this 
a measurement from an earlier monitoring campaign at Wise Observatory (Netzer et al. 1990), 
adjusted as described in Paper XV. The best-fit power-law relationship between the continuum 
flux and time lag is r cen t ot Fopt^i nearly a linear relationship. 

This result can be compared with a simple theoretical prediction. Photoionization equilibrium 
models for the BLR are characterized by (a) the shape of the ionizing continuum incident upon 
the BLR gas, and (b) an ionization parameter U, which is the ratio of ionizing photon density to 
particle density at the cloud inner face, 

U = 9(H) , (4) 

where (J(H) is the rate at which ionizing photons are produced by the continuum source, r is 
the separation between the BLR gas and the continuum source, and «h is the BLR particle 
density. Since the luminosity in ionizing photons L\ on is proportional to Q(K), we expect that 
r 2 oc Li on /Unn. Essentially, as Lj on varies, we expect that the response in a particular emission 
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line will be greatest at some particular value of the product Una, which will thus lead to r oc L ion 
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for any given line; in other words, the response of a particular line to a continuum variation ought 
to be dominated by gas in which the change in emissivity is greatest. In Figure 3, we show a test 
of this naive prediction based on fitting a r oc L 1 / 2 (i.e., r cent oc F^) power law to these data. 
The resulting fit is quite poor. 

The simple model above makes the implicit assumption that the continuum does not change 
shape as it varies. This allows us to use L opt as a surrogate for L lon , which cannot be measured 
directly. However, a flux at an ultraviolet wavelength closer to the Lyman edge at 912 A would 
obviously provide a better surrogate for L mn than L opt , and such data do exist in archival form, 
mostly from our previous UV monitoring experiments (Papers I and VIII) on NGC 5548. We have 
therefore recovered the NEWSIPS-extracted IUE spectra of NGC 5548 obtained since late 1989 
from the IUE Final Archive in order to effect a direct comparison of the UV and optical amplitude 
of continuum variability. We measured the UV continuum flux Fuv = -Fa (1350 A) by averaging 
the flux in the observed wavelength range 1370-1380 A for all IUE SWP spectra obtained within 
one day of optical continuum measurements. All optical fluxes within one day of a UV continuum 
measurement were used to form a weighted average optical flux. This procedure yielded 83 pairs 
of UV/optical continuum measurements. We then fitted a power-law function F opt oc Ffjy to these 
data, as shown in Figure 4, yielding a best-fit slope a = 0.56. Combining this with the above 
relationship between T cent and the optical continuum yields r ccn t oc Fyy 3 , which is consistent with 
the simple model. We note that there is no statistically significant time delay between the UV and 
optical continuum variations (Peterson et al. 1998b, and previous papers in this series). 

Thus, the apparently linear correlation between the variations in the optical continuum and 
those in the H/? emission line is attributable to similar relationships between the UV and optical 
continuum variations on the one hand and the UV continuum and H/3 emission-line variations 
on the other. This probably accounts for the apparent absence of scatter (relative to, say, Lya 
or Civ A1549) in comparisons of the optical continuum and H/3 fluxes of AGNs (e.g., Yee 1980; 
Peterson 1997). 

It is interesting that our observational result agrees with the naive theory once we account 
for the change in the shape of the continuum between the UV and optical, but that naive theory 
supposes that the shape of the ionizing continuum remains unchanged. Either the shape of the 
continuum between the UV and extreme ultraviolet still remains constant, or the H/3 line is 
surprisingly insensitive to the shape of the continuum. 

It should be pointed out that our result differs from the BLR radius-luminosity relationship 
found by Kaspi et al. (2000), namely r oc L 7 . The Kaspi et al. result describes how the BLR 
radius varies from object-to-object as a function of the mean optical luminosity of the source. The 
relationship discussed here describes how the BLR radius in an individual object changes as the 
luminosity of the central source varies with time. 
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5. Summary 

We have presented an additional five years (1997-2001) of optical spectroscopic and 
photometric observations of the continuum and H/3 emission-line variations in the Seyfert 1 
galaxy NGC 5548. We have also added significantly to the previously published data from 1996, 
warranting a complete recalibration of the data for that particular year. The new data expand 
our temporal coverage of variations in this AGN to a total of 1530 continuum and 1248 H/3 
measurements obtained over a 4756-day span, beginning in 1988 December and terminating in 
2001 December. During this period, the nonstellar optical continuum in NGC 5548 varied by 
approximately an order of magnitude in flux. 

Analysis of the time delay between continuum variations and H/3 response allows us to confirm 
our earlier finding (Paper XV) that the H/? lag varies with the mean continuum flux, ranging from 
a low value of ~ 6 days in 2000 to a high value of ~ 26 days in 1998-1999. We find consistency 
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with the simple photoionization equilibrium prediction r ce nt oc £ UV ' 
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Fig. 1. — The optical continuum (upper panel) and H/3 (lower panel) light curves from 1989 
December to 2001 December. The data are comprised of 1530 continuum measurements and 1248 
emission-line measurements. The continuum fluxes are in units of 10~ 15 ergs s _1 cm _2 A _1 , and 
the line fluxes are in units of 10~ 13 ergs s _1 cm -2 . The horizontal dashed line is an estimate of 
the continuum contribution from the host galaxy (Romanishin et al. 1995) through the standard 
aperture used here (5".0 x 7".5). Flux measurements are in the observer's reference frame and are 
uncorrected for Galactic extinction. 



Fig. 2. — Cross-correlation functions for the continuum and H/3 emission-line light curves of NGC 
5548. The various panels represent subsets of the new data presented here. The solid line shows 
the ICCF, and the vertical lines represent the la uncertainties associated with the DCF values. 
The peaks and centroids of the ICCFs are summarized in Table 8. 



Fig. 3. — The cross-correlation centroid r cen t as a function the yearly mean optical continuum flux 
(Fa (5100 A)), based on eight years of International AGN Watch data (filled circles) plus 1988 data 
from Wise Observatory (open circle). The vertical dashed line indicates our best estimate of the 
host galaxy flux F ga \ through our standard aperture. The solid line is the best-fit power law to the 
relationship r cen t oc F" t , where -F op t = -F\(5100 A) — F ga \, i.e., the nonstellar optical continuum 
flux. We find a best-fit slope a = 0.95. The dotted line is the best fit to the naive prediction 

Tcent OC -^opt ■ 



Fig. 4. — The relationship between UV continuum fluxes -Ftjv measured from NEWSIPS-extracted 
IUE spectra and the nonstellar optical continuum F opt for pairs of UV and optical continuum 
measurements made within one day of each other. Shown as a dashed line is the best-fit power 
law to these data, F opt oc F°f . Comb ining this relationship with that shown in Figure 4 yields 
Tccnt oc F^jy 3 , consistent with the simple model prediction T cen t oc F^Jy . 
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Table 1. Summary of Program on NGC 5548 



Year 


Dates 




Comments 


(!) 


(2) 




( 3 ) 


1 


1988 Der 14 - 1989 Aue 7 


1; 


TTV pontmnnTYi ann lmp^ 

\J V L L*ll L 111 LL Lllll CI 1 ILL 1111CO 




1988 Dpp 14 - 1989 Opt 10 


9- 


( ) t ~> f" i /->Q 1 pnntinn 1 1 tyi q n rl TT f{ 

KJ jj L1L,CX1 L,UllLlllLi Lllll CXllLL L±j^J 






s- 

KJ • 


A/Tnlti'nlp OTvtipal linpt; 

1\± 111 LI JJ1L- VJ|JL1L/Cll 1111CO 






4- 


T^almPT* Cf~\Y\ ll Till 1 1 TY1 Q n H hIp TT 
J_JdllllCl LvUllljlllLi LL111 ClllLL ± v_, 11 






5- 


( ity!~i pa 1 i tyi a on tics' / Ti n ntnTTi p^tv 

V-' L* LlL/CLl llllCltilllci / U11LJ ULJ111C Li y 






6- 


T? pvi tipn r»"nT"i pa 1 pont i nn 1 1 tyi 

IICVIOCU VJ 1J L1L,C%1 L-LJ11L111LL Lllll 


2 


1989 Dpp 1 - 1990 Opt 15 


fi- 


( )"nt i pa 1 pont i nn 1 1 tyi a n rl TT H 

KJ b LVjCXil L-LJ11L111LL Lllll ClllLl xLj-.J 




1 990 Nov 29 - 1 991 Opt 5 


7- 


OirHpal pnntimniTYi anrl TT/? 

KJ JJUlL,Clil LvLJllLlllLL Lllll ClllLL ±L^J 


4 


1992 Tan 1 - 1992 Opt 3 


7- 


OtyHppiI pontinninn anrl TT/? 

V-/ LlL-CLl L-LJ11L111LL Lllll ClllLl ±±1^/ 


5 


1993 Mar 14 - 1993 Mav 27 


8: 


TIV pont.innum and lines 

KJ V ^ V_7 llulll LL l_l All CAill V_l llll^O 




1992 Nov 27 - 1993 Sep 25 


8: 


Optical coiitinuuin and H/? 




1993 Mar 10 - 1993 May 14 


9- 


HIT T \/ rnntinn 1 1 vn 

1— - L. 1 V L/UilUlllLL Lllll 


6 


1993 Nov 17 - 1994 Oct 11 


10 


■ Ontical continuum arid TT/? 

t \ — r l_y LflV^LLil VjVJllUlllLi Lllll UillLl -L-LA^/ 


7 


1994 Nov 22 - 1995 Oct 17 


10 


* Ontical continuum and TT/? 

i ^ — f Y-J VI ± V.' Ll!l V_/W11U111L1 Lllll CLillLA -L-LA^/ 


8 


1995 Nov 22 - 1996 Oct 16 


10 


: Optical coiitinuuin and H/3 






11 


■ T? pvi Qpn nntiral rnntirin 1 1 tyi ann TT H 

1LC V1DCU UJJ LlL-CLl L-U11L111LL Lllll ClllLL ±Lj^J 


9 


1996 Dec 17 - 1997 Oct 7 


11 


( )ty1~i pa 1 print inn 1 1 tyi a n n TT fi 

• V^/ VJ L1L-CL1 L-VJllUlllLl Lllll ClllLL J. l-kJ 


10 


1997 Nov 22 - 1998 Sep 28 


11 


\ Optical continuum and H/3 


11 


1998 Nov 24 - 1999 Oct 4 


11 


■ Ontical continuum and TT/? 

» V ^ l_/ LflV^Llil VjL/11L111LI Lllll CL11L1 Illy 


12 


1999 Dec 5 - 2000 Sep 3 


11 


■ Ontical continuum and TT/? 

* ^ — ^ l_/ LflV^LLil VjV71HJ111LI Lllll UillLl J-J-f^s 


13 


2000 Nov 30 - 2001 Dec 21 


11 


■ OirHpal pnntinnnm anrl TT/? 

KJ JJ U1L-CL1 L/VJ11U111L1 Lllll ClllLL ±±j^S 


References: 






1 


Clavel et al. 1991 (Paper I) 






2 


Peterson et al. 1991 (Paper II) 






3 


Dietrich et al. 1993. (Paper IV) 






4 


Maoz et al. 1993 






5 


Romanishin et al. 1995 






6 


Peterson et al. 1992 (Paper III) 






7 


Peterson et al. 1994 (Paper VII) 






8 


Korista et al. 1995 (Paper VIII) 






9 


Marshall et al. 1997 






10 


Peterson et al. 1999 (Paper XV) 







11 This work. 
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Table 2. Sources of Spectroscopic Observations, Years 8-13 


Data 




Number of 


Set 


Telescope and Instrument 


Spectra 


(1) 


(2) 


(3) 


A 


1.8-m Perkins Telescope + Ohio State CCD Spectrograph 


57 


F 


1.5-m Mt. Hopkins Telescope + CCD Spectrograph 


258 


GH 


2.1-m Guillermo Haro Observatory Telescope + B&C Spectrograph 


22 


H 


3.0-m Lick Shane Telescope + Kast Spectrograph 


28 


K 


2.4-m MDM Observatory Hiltner Telescope + MODSPEC Spectrograph 


6 


L 


6-m Special Astrophys. Obs. Telescope + CCD Spectrographs 


19 


LI 


1-m Special Astrophys. Obs. Telescope + CCD Spectrographs 


32 


M 


3.5-m and 2.2-m Calar Alto Telescopes + CCD Spectrograph 


7 


R 


1.5-m Loiano Telescope + CCD Spectrograph 


1 


W 


2. 6-m Shajn Telescope + CCD Spectrograph 


113 




TOTAL 


543 



Table 3. Flux Scale Factors for Optical Spectra 



Data Point-Source Extended Source 

Set Scale Factor Correction G 

cp (1CT 15 ergs s _1 cm~ 2 A -1 ) 
(1) (2) (3) 



A 

F (Year 8) 

F (Years 9-13) 

GH 

H 

K 

L 

LI (Year 10) 
LI (Year 11) 
LI (Year 12) 
LI (Year 13) 
M 
R 

W (Year 8) 
W (Years 9-13) 



1 


.000 











.941 


± 


0. 


.035 


1 


.039 


± 


0. 


,042 





.985 


± 





.032 


1 


.013 


± 


0. 


,044 


1 


.010 


± 


0. 


,050 


1 


.003 


± 


0. 


,056 


1 


.113 


± 


0. 


.045 


1 


.011 


± 


0. 


,024 


1 


.054 


± 


0. 


,064 


1 


.011 


± 


0. 


.094 





.990 


± 





.014 





.863 


± 





.025 





.932 


± 





.031 


1 


.009 


± 


0. 


.022 



0.000 

-1.730 ±0.050 
-0.532 ±0.570 
-0.329 ±0.555 
-0.418 ±0.404 
0.234 ± 0.448 
-0.874 ± 0.645 
4.088 ± 0.396 
1.500 ±0.931 
2.556 ± 0.844 
3.090 ± 0.954 
-1.353 ±0.380 
-1.900 ± 1.065 
-1.319 ±0.437 
-0.422 ± 0.244 
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Table 4. Sources of Photometric Observations, Years 9-13 



Data 




Number of 


Set 


Telescope and Instrument 


Observations 


(1) 


(2) 


(3) 


A 


0.6-m Telescope, Crimean Laboratory of the Sternberg 
Astronomical Institute + pulse-counting 
UBV photometer, 14".3 aperture. 


91 


B 


1.25-m Telescope, Crimean Astrophysical Observatory 
+ pulse-counting Photometer-Polarimeter, 15" aperture 


12 


C 


1.56-m Telescope, Shanghai Astronomical Observatory 
+ CCD, 10" aperture. 


7 


D 


0.48-m and 0.6m Telescopes of Ulugbek Astronomical 
Institute + pulse-counting UBV photometer, 14".3 aperture. 


69 


E 


1-m and 0.6-m Telescopes of the Special Astrophysical 
Observatory + CCD, 10" aperture. 


67 



TOTAL 246 
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Table 5. Optical Continuum and H/3 Light Curves' 



Julian Date 


F A (5100 A) 


F(EP) 


(-2400000) 


(10~ 15 ergs s _1 cm" 2 A _1 ) 


(10~ 13 ergs s _1 cm -2 ) 


(1) 


(2) 


(3) 


50044.05 


11.52 ±0.23 


9.56 ±0.19 


50048.60 


10.97 ±0.49 


8.78 ±0.23 


50052.08 


10.77 ±0.27 


8.96 ±0.31 


50053.01 


10.33 ±0.21 


9.02 ±0.18 


50061.02 


8.56 ±0.17 


9.19 ±0.18 


50064.60 


8.20 ±0.37 


8.33 ±0.22 


50069.09 


8.18 ±0.20 


8.22 ±0.29 


50074.00 


8.18 ±0.33 


7.27 ±0.20 


50075.00 


8.30 ±0.33 


7.04 ±0.20 


50078.00 


8.15 ±0.33 


6.92 ±0.19 



a A complete version of Table 5 is available at 



http: / /www. astronomy.ohio-state.edu/ ^agnwatch / 
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Table 6. 


Sampling Statistics for Optical Continuum 








Number of 


Sampling 


Interval (days) 


Mean 






Subset 


Epochs 


Average 


Median 


Flux a 


- 1 var 


-Rmax 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


All data 


1530 


3.1 


1.0 


9.73 ± 2.44 


0.247 


3.57 ±0.18 


Year 1 (1989) 


125 


2.4 


1.0 


9.92 ± 1.26 


0.117 


2.16 ±0.16 


Year 2 (1990) 


94 


3.4 


2.0 


7.25 ± 1.00 


0.129 


1.82 ±0.09 


Year 3 (1991) 


65 


4.8 


3.0 


9.40 ± 0.93 


0.090 


1.51 ±0.09 


Year 4 (1992) 


83 


3.4 


2.0 


6.72 ± 1.17 


0.168 


2.04 ±0.10 


Year 5 (1993) 


174 


1.3 


0.7 


9.04 ±0.90 


0.092 


1.65 ±0.08 


Year 6 (1994) 


135 


2.4 


1.0 


9.76 ± 1.10 


0.104 


1.76 ±0.12 


Year 7 (1995) 


83 


4.0 


1.9 


12.09 ± 1.00 


0.079 


1.48 ± 0.04 


Year 8 (1996 


144 


2.3 


1.0 


10.56 ± 1.64 


0.150 


1.86 ±0.11 


Year 9 (1997) 


126 


2.4 


1.0 


8.12 ±0.91 


0.105 


1.84 ±0.06 


Year 10 (1998) 


175 


1.8 


0.6 


13.47 ± 1.45 


0.100 


1.59 ±0.09 


Year 11 (1999) 


148 


2.1 


1.0 


11.83 ± 1.82 


0.149 


1.98 ±0.11 


Year 12 (2000) 


94 


2.9 


1.0 


6.98 ± 1.20 


0.166 


2.40 ±0.12 


Year 13 (2001) 


84 


4.6 


2.0 


7.03 ± 0.86 


0.112 


1.68 ±0.10 



a Units of 10 15 ergs s 1 cm 2 A 1 . 
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Table 7. 


Sampling 


Statistics for H/3 Emission Line 








Number of 


Samplin 


g Interval (days) 


Mean 






Subset 


Epochs 


Average 


Median 


Flux a 


L var 


^max 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


All data 


1248 


3.7 


1.1 


7.82 ± 1.71 


0.216 


4.60 ± 0.34 


Year 1 (1989) 


132 


2.3 


1.0 


8.62 ±0.85 


0.091 


1.57 ±0.12 


Year 2 (1990) 


94 


3.4 


2.0 


5.98 ± 1.17 


0.191 


2.30 ±0.12 


Year 3 (1991) 


65 


4.8 


3.0 


7.46 ± 0.81 


0.093 


1.58 ±0.14 


Year 4 (1992) 


83 


3.4 


2.0 


4.96 ± 1.44 


0.284 


3.03 ± 0.30 


Year 5 (1993) 


142 


2.1 


1.0 


7.93 ± 0.53 


0.057 


1.40 ± 0.06 


Year 6 (1994) 


128 


2.6 


1.0 


7.58 ± 0.94 


0.117 


1.57 ±0.07 


Year 7 (1995) 


78 


4.2 


2.1 


9.27 ±0.70 


0.071 


1.38 ±0.06 


Year 8 (1996) 


144 


2.3 


1.0 


7.95 ± 0.87 


0.106 


1.76 ± 0.08 


Year 9 (1997) 


95 


3.1 


1.1 


7.41 ± 0.95 


0.125 


1.66 ±0.07 


Year 10 (1998) 


119 


2.6 


1.0 


10.27 ± 1.03 


0.097 


1.50 ±0.07 


Year 11 (1999) 


86 


3.7 


1.0 


9.34 ±0.61 


0.058 


1.50 ±0.06 


Year 12 (2000) 


37 


7.6 


3.0 


6.27 ± 1.22 


0.192 


1.96 ±0.07 


Year 13 (2001) 


45 


6.7 


3.2 


5.26 ± 1.12 


0.208 


2.33 ± 0.08 



a Units of 10 13 ergs s 1 cm 2 . 
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Table 8. Cross-Correlation Results 







''"cent 


Tpeak 




Subset 


(days) 


(days) 


''max 




(1) 


(2) 


(3) 


(4) 


Year 


1 (1989) 


19.73l 2 -°J 


21.511-? 


0.877 


Year 


2 (1990) 


19.34+H 


i8.5^:? 


0.909 


Year 


3 (1991) 


16.35_ 3 28 


17.5l|? 


0.740 


Year 


4 (1992) 


1 1 07+2.30 
-2.30 


13.712"? 


0.918 


Year 


5 (1993) 


13.61±J;£ 


13.2l 3 | 


0.730 


Year 


6 (1994) 


1 r 4q +2.26 


8.7l|| 


0.832 


Year 


7 (1995) 


21-43t|^ 


22.911} 


0.880 


Year 


8 (1996) 




15.31^ 


0.907 


Year 


9 (1997) 


17.49l 2 /^ 


17 7 +3.4 
1 '-'-1.8 


0.803 


Year 


10 (1998) 


26.44±|| 


01 7 +0.5 
oi -'-10.9 


0.844 


Year 


11 (1999) 


25.75_ 3 ; 32 


oa 4+1-0 
— 18.5 


0.846 


Year 


12 (2000) 


D - iz -4.08 




0.874 


Year 


13 (2001) 


14.89111? 


8.0li 4 4 6 


0.793 




l i i i l i i i i l i i i i l i i i i l i i i i l i 
48000 49000 50000 51000 52000 

JD - 2,400,000 




Lag (days) 




Mean Flux F x (5100 A) (10 15 ergs s" 1 cm" 2 A" 1 ) 



